Spectral analysis of the gamma-ray background near the dwarf Milky Way satellite 
Segue 1: Improved limits on the cross section of neutralino dark matter annihilation 
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The indirect detection of dark matter requires that dark matter annihilation products be discrim- 
inated from conventional astrophysical backgrounds. Here, we re-analyze GeV-band gamma-ray 
observations of the prominent Milky Way dwarf satellite galaxy Segue 1, for which the expected 
astrophysical background is minimal. We explicitly account for the angular extent of the conserva- 
tively expected gamma-ray signal and keep the uncertainty in the dark-matter profile external to 
the likelihood analysis of the gamma-ray data. 

PACS numbers: 95.35.-|-d; 95.85.Pw; 95.30.Cq 



I. INTRODUCTION 

The physical nature of dark matter is still mysteri- 
ous. One of the most probable hypotheses is that it is a 
mixture of equal quantities of weakly-interacting massive 
particles (WIMPs) and their antiparticles, generated in 
the early Universe. The detection of their annihilation 
products could be an attractive possibility to clarify the 
nature of these particles. 

Dwarf satellites galaxies of the Milky Way are a class 
of objects that are very promising for detecting dark- 
matter annihilation. Though the dark-matter density 
in these objects is significantly lower than that at the 
Galactic Center, they typically contain no sources of cos- 
mic rays, and therefore almost any high-energy signal 
detected from the dwarf galaxies could be the signal of 
dark-matter annihilation. 

Tens of Milky- Way satellites are known to date. In 
our work we examine the dwarf galaxy Segue 1, since 
it combines several advantages: it is the closest satellite 
to Earth, it is located at high Galactic latitude (and, 
consequently, low gamma-ray background), and among 
the dwarf galaxies. Segue 1 is expected to produce the 
strongest dark- matter signal [ll-Il]. Due to it's proximity 
to the Sagittarius stream, the nature of Segue 1 has been 
disputed: it was argued to be a disrupted star cluster 
originally associated with the Sagittarius dSph [3| . How- 
ever, a kinematic study of a larger member-star sample 
(66 stars compared to the previous 24-star sample) has 
recently confirmed that Segue 1 is indeed a dwarf satellite 
galaxy Q. Table 1 in Q shows that the line-of-sight an- 
nihilation integral through the dark matter distribution 
of Segue 1 integrated over the angular size of the source 
is the largest of all Milky Way dwarf satellites. 

Other approaches used in the literature include per- 
forming a joint likelihood analysis of several dwarf galax- 
ies and a Bayesian technique 0. We do not find 
a joint analysis advantageous: actually, only a few ob- 
jects can produce a flux comparable to Segue 1. There- 
fore, the statistics can only be improved marginally, 
while the merging of several sources with various prop- 



erties, density-profile measurement errors, and back- 
grounds strongly increases the ambiguities of the anal- 
ysis. 

Our consideration is based on the assumption that the 
WIMP in question is the lightest SUSY particle, a neu- 
tralino. However, the only property of neutralinos we use 
is that they annihilate mainly into a pair of heavy par- 
ticles (heavy quarks, gauge bosons, tau leptons etc.). So 
our conclusions are valid for any WIMPs for which the 
supposition applies. 

We use data from the Fermi Large Area Telescope 
(LAT), since this telescope is currently the most sensi- 
tive to the photons that can be generated by the annihi- 
lation of neutralinos with mass roughly between 10 GcV 
and 1 TeV. Currently operating imaging Cherenkov tele- 
scopes are primarily sensitive to gamma rays with ener- 
gies > 200 GeV and hence constrain very massive dark- 
matter particles [1^ . Jermi-LAT data have already been 
used in several articles 0, [1,0, Hi- However, all of them, 
except for Q, considered Segue 1 as a point source. We 
perform an accurate spectral analysis and explicitly ac- 
count for the spatial extent of the source, using more data 
than 0, which allows us to improve the result. In jl|, 
the Segue- 1 dark- matter distribution is modeled with a 
Navarro-Frenk- White profile 9| , a centrally-concentrated 
density distribution, thus leading to an enhanced gamma- 
ray signal. We use a template for the angular distribution 
based on the significantly more conservative Einasto pro- 
file with arbitrary scale factor, diffuse and point-source 
backgrounds. The Einasto profile predicts a spatially 
more extended signal, and the nominal annihilation flux 
is weaker than in Q . As we will show in section |Vl 
the Einasto model predicts an annihilation signal that 
is only marginally higher than the lowest limit for any 
density distribution compatible with the stellar dynam- 
ics of Segue 1. We further investigate a possible position- 
dependent Sommerfeld enhancement [lO| of the annihi- 
lation rate. As both the total annihilation signal and 
the source extent are affected by the choice of dark- 
matter profile and enhancement factors, the uncertain- 
ties project onto the gamma-ray based upper limits in a 
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FIG. 1: The angular distribution of the signal for the four 
profiles used in this article. Here 9 is the angular distance 
from the center of the source, F{r) is defined by ([2} using 
C = 1 and {(jv) — {av)c- 

non-trivial way. 

II. MODEL DESCRIPTION 

In order to avoid ambiguity, we hereafter suppose that 
the dark-matter particle is identical to the antiparticle. 
If the WIMP is a Dirac particle (which is not the case 
for neutralinos), then the annihilation rate halves and all 
limits on the cross-sections would be weaker by a factor 
2. On the other hand, the main purpose of our work is to 
compare the limits obtained with the cross-section esti- 
mate from dark-matter abundance {(Jv)a — 3 • 10~^^ cm'^ 
s"""^ [m (see below). This value is also increased by the 
same factor 2, if the dark-matter particle is not a Ma- 
jorana particle. Therefore, the ratio between the limit 
obtained and {av)a is not sensitive to the nature of the 
dark matter in this sense. 

The annihilation signal (i.e., the number of photons d'^yi 
arriving at the local observer from a solid angle element 
dil, per unit time interval dt, per an area dA, and an 
energy interval dE) in the general case is 

^_ dm _ jav), ^ dN 

dAdfldEdt Snml dE ^ ' 

where dN/dE is the spectrum of the photons produced 
by a single annihilation event, and 

F{r,) = / C^pHl (2) 

JlOS {'^V)c 

F{rp) is the line-of-sight integral of the multiplication of 
tree factors through the density profile of the source at 
the projected radius Vp. C = (p^)/(/o)^ is the so-called 
boost factor, which takes into account the enhancement 
of the annihilation signal by dark-matter substructures 



that may be present. The question of possible substruc- 
tures in Segue 1 is extremely vague, and we will not con- 
sider it, accepting C = 1. The factor {av)/{av)c al- 
lows for a possible dependence of {av) on particle speed, 
e.g. the Sommerfeld effect. {(Jv)c is a characteristic 
value chosen such that for absent velocity dependence 
{av) ~ {iyv)c — constant. 

With knowledge of the distance, F{rp) fully describes 
the angular distribution. We should mention that quan- 
tities such as J(n) = L{n) = F {vp [fi/]) d^l/, where 
J or L are considered as a function of the angle, are of 
more common use in the literature. Being defined so, 
they depend on the distance to the source, contrary to 
F{rp). We set the distance to Segue 1 to 23 kpc. 

The mass of Segue 1 within a certain radius is mea- 
sured by observations of a few (< 100) member stars, 
but the density profile is assumed [e.g. Q, [T^- To de- 
scribe the signal profile we used four models. The first 
two assume that (av) does not depend on the velocity. 

The first model (hereafter the standard model) is based 
on the Einasto density profile: 
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with the parameters of Equation|3]in accordance with |13| 
as ps = 1.1 • 10® Mokpc"^, rs = 0.15 kpc, and n = 3.3. 
Further we assume for the standard profile (av) — {av)c. 
It is worth mentioning that the enclosed mass and the 
function F{rp) can be experimentally defined much bet- 
ter than the individual parameters p^, Tj, and n. Some 
authors |7| use a completely different set of parameters, 
which corresponds to a quite similar F(rp). We note that 
the total integral Jtot = § Fdfl = lO^^ ^ GeV^srcm-^ 
in our case (and cgs units) is smaller than the Jtot — 
10^^-^ GeV^srcm-s used in [|. 

We should emphasize significant uncertainties of the 
Segue 1 density profile: almost all the stars lie inside 
10' (~ 67 pc) from the center [3], which is significantly 
smaller than rg. In order to test how sensitive the an- 
nihilation signal is to the choice of profile, we used a 
second (minimal) model. The assumptions of the min- 
imal model are that the density distribution beyond a 
radius of 67 pc and the total mass inside this radius co- 
incide with those of the standard model ([3]), but the dark 
matter is homogeneously distributed in the inner 67 pc. 
In this case the predicted signal is minimal among the 
profiles compatible with observations. This conservative 
choice of density profile also renders the angular extent of 
the annihilation signal larger than in the case of a more 
spiked density profile, for which a point-source analysis 
may be suitable. 

Among all the scenarios with (av ) depending on v, we 
consider the Sommerfeld effect [l3|- If the dark-matter 
particles interact only via standard gauge bosons, the ef- 
fect arises only for heavy WIMPs, with masses in the 
TeV range, but if there are lighter mediator particles, 
the enhancement appears even in the case of low-mass 
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FIG. 2: The photon spectra generated by the annihilation of 
a pair of 200- GeV neutrahnos through various channels (the 
spectra are normalized as if the channel under consideration 
were the only channel of annihilation). The bb, cc, rf, tt, 
WW, and ZZ are represented by the lines with squares, cir- 
cles, stars, triangles, W, and Z, respectively. 

WIMPs • A similar explanation has been invoked for 
the PAMELA positron excess however, any discus- 
sion of the PAMELA results or of the boosting necessary 
for the dark matter interpretation of the PAMELA ex- 
cess goes far beyond the scope of this article. We would 
like to emphasize only two essential points. First, the 
Sommerfeld effect can take place only if there is some 
new relatively light mediator field for the majority of the 
WIMP masses under consideration in this article (except 
for the highest ones). Second, if the Sommerfeld effect 
takes place, (av) depends on the collision velocity, and 
we should indicate the characteristic velocity when we 
obtain experimental constrains of the cross-section. 

To account for the Sommerfeld effect we used two mod- 
els. Both of them are based on the Einasto density pro- 
file, as the standard model. This allows us to find (f>{r), 
the gravitational potential of Segue 1. Roughly speaking, 
if the Sommerfeld effect occurs, then (av) ^ 1/v. The 
question is how v depends on radius, r. We scale (av) 
using as characteristic velocity the escape velocity from 
the center of Segue 1, Wc = 2 ^—0(0) ~ 60 km s^^: 

{av) = ^{av),. (4) 

V 

For both the models taking into account the Sommerfeld 
effect the limits on the cross-section we obtain are con- 
sidered to mean the constraints on the quantity {av)c, 
the {(Tv) at the average collision velocity Vc- {<Jv) at the 
decoupling epoch could be much lower in this case. 

Besides the absolute intensity, the Sommerfeld effect 
affects the angular profile: we should take into account 
the velocity distribution of the particles as a function of 
r. The velocity distribution properties are still not quite 
clear [l^. We consider two possible hypotheses about 



the distribution of average collision speed v inside Segue 
I. We may suppose that the speed is proportional to 
the escape velocity, or, for simplicity, v — 2vesc/V^ = 
2\J—(j){r). A very similar situation takes place if radial 
motion of the particles dominates ^7| . Hereafter we call 
this model virial. 

In the opposite case, when the distribution function is 
supposed to be more or less isotropic, it is reasonable to 
suppose V = ^/2vorb{r), where Vorbir) is the orbital speed 
at radius r. This situation resembles the well-known 
isothermal dark-matter halo and should be typical for 
the halo center [l^. Moreover, the experimentally ob- 
served small velocity dispersion of stars at the center 
of Segue 1, ~ 3.7 km s~^, favors this scenario. We name 
this model orbital. 

The signal profile as described by each of the four mod- 
els is shown in Figure [TJ 

The error in the experimental determination of the 
multiplier F(rp) slightly depends on r and can be es- 
timated [3| as 

a(logioF)~0.6 . (5) 

If we assume the uncertainty in logj^g ^o be Gaussian, 
then the 2-a confidence interval of Fij-p) corresponds to 
an additional multiplier in equation ([T]) extending over 
the range 0.063 — 16. It is obvious that this is the main 
source of uncertainty in the determination of the upper 
limit on the cross section. However, F(rp) is an external 
factor in the likelihood analysis of the fermi-LAT data, 
provided the angular profile is realistic, and therefore for 
a given dark-matter profile and signal boosting the limits 
on {av) will to first order scale with the error factor in 

nrp)- 

The angular size of the source is comparable to the 
angular resolution of Fermi-LAT above 10 GeV; for the 
standard model 68% of the signal is expected within 0.17° 
and 95% of the signal within 0.54° of the source cen- 
ter. As we used a binned likelihood analysis, we wanted 
to avoid a dependence of the model on the position in 
the map grid. For that purpose we smoothed the pro- 
files: Instead of F{rp) defined by ^ we used F{rp) = 
/ao F{rp)dn/Afl, where AQ is ~ 0.037° x 0.037°, which 
is slightly smaller than the pixel size and corresponds to 
an area 15pc x 15pc at the source. The smoothing kernel 
is about a factor 2.5 narrower than the angular resolution 
of Fermi-LAT above 10 GeV and its particular choice has 
no impact on the results of the likelihood fit. 

The question of the photon spectrum produced by neu- 
tralino annihilation is, generally speaking, quite complex. 
Annihilation into a pair of particles (quarks, leptons, 
gauge bosons etc.) is the most likely process. The chan- 
nels with light final states are suppressed by helicity, and 
the primary products of the process are heavy particles 
(6, t, c-quarks, W, Z-bosons, r leptons). Let us denote 
the branching ratio of each individual channel by k (for 
instance, krr is the average fraction of annihilations re- 
sulting in a TT pair generation). The heavy primary par- 
ticles are unstable and decay into stable products includ- 
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ing photons (mainly from pion decays). We calculated 
the spectrum of photons generated by various annihila- 
tion channels with the DarkSUSY package [l9l |. 

The real photon spectrum is a linear combination of 
the spectra corresponding to separate channels with coef- 
ficients fcj. Channel contributions ki (and, consequently, 
the spectrum) depend on the choice of SUSY model. To 
avoid a complete scan of the SUSY parameter space, we 
use the spectra of separate channels to obtain upper lim- 
its on the products fc^cr. For each neutralino mass we use 
two types of spectral models: the spectra corresponding 
to pure annihilation into bb and rf, respectively. We re- 
produce in Fig. [2] the photon spectra generated by the 
annihilation of a pair of 200-GeV neutralinos through 
various channels (the spectra are normalized as if the 
channel under consideration were the only channel of an- 
nihilation). As we can see, the spectra of all channels, 
except for rf and tt, are very similar to that of the bb 
channel. Therefore, by placing a limit on khtcr, we ac- 
tually provide a very good approximation of the upper 
bound on {kbb + kcc + kww + kzz) o". 

Note that, for a significant part of SUSY parameter 
space, (fcfib + kcc + kww + kzz) is close to 1. The relative 
width of the rf channel is usually significantly smaller 
than 1, but the spectrum is very hard and, consequently, 
easily distinguished, since the background rapidly de- 
creases with energy. The tt spectrum is somewhat softer 
than that of bb decay, but we choose not to consider 
it as the third spectral model for a number of reasons. 
First of all, this channel is closed for light neutrali- 
nos because of the large mass of the i-quark. Second, 
(fcbf, -I- kcn + kww + kzz) is typically at least comparable 
with ktt [20] . Third, the spectrum of the tt channel does 
not differ very much from that of bb. 

We consider an array of neutralino masses that covers 
the expected mass range [21[ approximately uniformly 
in \ogE: 40, 50, 63, 80, 100, 126, 159, 199, 251, 316, 
398, 501, 631, 794, 1000, 1259, 1585, 1995, 2512, 3162, 
3981, and 5011 GeV. In total we use 44 spectral models: 
a hh and a rf model for each of 22 neutralino masses. 
The spatial models are, of course, independent on the 
neutralino mass. 



III. OBSERVATIONS OF SEGUE 1 WITH 
FERMI-I.AT 

We searched for evidence of gamma-ray emission from 
the dwarf spheroidal galaxy Segue 1 using data from 
LAT, the principal instrument on board the Fermi space- 
craft, a description of which is given in [22]. 

Data analysis was performed with the Fermi-LAT 
ScienceTools version 9.23.1, using the post-launch 
P7SOURCE_V6 instrument response functions. The 
cosmic-ray background was reduced by choosing the 
P7S0URCE class which only includes events with high 
probability of being photons. We also require a zenith- 
angle cut of 100° to minimize the spill-over of photons 
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FIG. 3: 95%-C.L. integral flux upper limits for Segue 1. The 
curves with the same color represent the two channels 100% 
bb (triangular markers) and 100% rf (square markers) for a 
specific model. All of the four models described in the text 
are shown. 



from the Earth limb, excluding time intervals where any 
part of the region of interest (ROI) was beyond the 
zenith-angle limit. In addition, data were not taken 
while the observatory was transiting the South Atlantic 
Anomaly and excluded when the rocking angle exceeded 
52°. The data used in this work came from observations 
in the period 2008 August 8 - 2012 February 4. 

A binned likelihood analysis is performed on a map 
with 0.05° pixel size in gnomonic (TAN) projection, mea- 
suring 10° on a side and centered on Q!J2ooo = 151.767° 
and 5j2ooo = 16.0819°, the nominal position of Segue 1. 
We only used photons with reconstructed energy greater 
than 2.4 GeV, for which the 68%-containment radius 
of the point-spread function (PSF) is narrower than 
^ 0.45°. The annihilation spectra are always substan- 
tially harder than those of the background, and therefore 
a high energy threshold can help to minimize background 
contamination. 

The resulting background model includes seven sources 
listed in the second LAT source catalog and 
the LAT standard Galactic diffuse emission component 
gal_2yearp7v6_v0 . f its along with the corresponding 
isotropic template iso_p7v6source.txt that accounts 
for extragalactic emission and residual cosmic-ray con- 
tamination [2^. All seven sources lie outside our ROI 
and are modelled as point sources with power-law energy 
spectra. Their spatial and spectral parameters were kept 
fixed at the values given in the catalog, while we per- 
mitted the normalizations of the diffuse components to 
freely vary. 

Segue 1 is characterized by a spatial template of 2° 
diameter with a resolution of 0.05°, following Equations 
[2] and [3l Despite of the fact that the spatial binning is 
slightly larger than the smoothing kernel applied to the 
Segue 1 profile, the results of the likelihood analysis are 
not affected. Analyses with different binning scales show. 
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in fact, changes in the integral flux upper limits less than 
1%. 

As described above, we use two different model spec- 
tra for the self-annihilation of neutralinos into quarks (bb) 
and leptons (rf ) for each of the four model profiles (stan- 
dard, minimal, virial, orbital). We tested 22 values of 
dark-matter particle mass in the range from 40 GeV to 
5011 GeV. 



m^, GcV 


Flux U.L., bb 


Flux U.L., Tf 


50 


7.98 


7.65 


100 


8.43 


5.83 


501 


7.42 


3.44 


1000 


6.64 


2.99 


5011 


5.24 


2.80 



TABLE I: Integral-flux upper limits above an energy thresh- 
old of 2.4 GeV for different WIMP masses (standard model). 
The confidence level is 95% and the units 10^^^ ph cm^^ s^^. 



IV. ANALYSIS AND RESULTS 

We analyzed data between 2.4 GeV and 300 GeV. The 
low-energy limit was a priori chosen to optimize the con- 
straints on the neutralino masses above 100 GeV. The 
normalization of the spatial template of Segue 1 is the 
relevant free parameter since we keep the position fixed. 
No significant gamma-ray signal was detected. We thus 
derive integral-fiux upper limits over the energy range 
2.4 GeV - 300 GeV using the profile likelihood tech- 
nique. Table H] provides 95% confidence-level upper limits 
for dark-matter self-annihilation emission from Segue 1 
(standard model only). 

Figure [3] shows the upper limits on the integrated fiux 
as a function of the neutralino mass for the two annihi- 
lation channels for the four models. We note that the 
bb final state predicts on average a smaller photon en- 
ergy than the rf final state (cf. Figure [2]). Therefore, 
the background contamination is more serious and, con- 
sequently, the fiux upper limits are lower for the rf chan- 
nel. 

To be noted from Figure 1 is the effect of the shape 
of the annihilation profile, F(rp). The 30-% difference 
between the best and worst upper limit reflect the dete- 
rioration in sensitivity arising from source extent. The 
flux upper limits do not depend on the normalization of 
the annihilation proflle, the astrophysical factor J. 

The impact of using P7S0URCE class events, in- 
stead of the older P6_V3_DIFFUSE, was also investi- 
gated. Differences of the two data sets slightly affect 
the results of the upper limits on the velocity-averaged 
self-annihilation cross section of neutralinos. The anal- 
ysis performed with the Pass-6 data provides more con- 
straining upper limits compared with the Pass-7 data, of 
around 10% for the low-mass neutralino (i.e., m^^, < 100 
GeV for bb and m^ < 501 GeV for rf model) and around 
20% for those with higher mass. Note that this change 
is likely due to a statistical fluctuation. The templates 
for diffuse and isotropic emission are different for Pass 6 
and Pass 7, and they dominate the likelihood fit. In the 
Pass-7 list of events between 2.4 GeV and 300 GeV re- 
constructed energy, we find 7 events within 0.5° of the 
nominal position of Segue 1 (5 events in the Pass-6 event 
list). The entire region of interest comprises 920 events, 
or 9.2 per square-degree on average. 

Comparing the standard profile with the point-source 
limit, we conclude that accounting for the source extent 
in the likelihood analysis increases the upper limits by 



typically 20%-25%, but only 10%-15% for rf models 
with large particle mass, for which most of the gamma 
rays have energies > 100 GeV. 



V. DISCUSSION 

Figure [5] represents the restrictions on (cru) imposed 
by our analysis, assuming kht = 1 or k-rr — 1- Note that 
the restrictions do not take into account the uncertainty 
(Eq. [5|) in the dark-matter distribution. The systematic 
and statistical uncertainties arising from modelling the 
stellar kinematics of Segue 1 are in fact the dominant 
source of error, for which we do not know the probability 
distribution. In any case, it is not a purely statistical 
uncertainty. 

We can draw several conclusions. First of all, the limits 
for the standard and minimal models differ only weakly, 
suggesting that the limits obtained using the Einasto pro- 
file are quite conservative and not very sensitive to the 
shape of the density profile. The limits for the standard 
and virial models almost coincide. This is a consequence 
of the fact that the central region of the source is rather 
small, and its angular size is smaller than r^. As a re- 
sult, (j> (and so the average particle speed) changes only 
slightly inside the source. 

The only model that excludes the relic abundance 
cross-section {av)a — 3 • 10^^^ cm^ (in the area of 
small neutralino masses) is the orbital. This is due to the 
fact that Vorb — >■ if the central density of the halo is 
finite. As a result, the signal in the center is significantly 
boosted. Unfortunately, the limitation is significant only 
for low-mass neutralinos and in the presence of the Som- 
merfeld enhancement. As we have already discussed the 
Sommerfeld effect for low-mass particles can appear only 
if the WIMPs interact via some low-mass mediator bo- 
son, discriminated from the known gauge bosons. So we 
can exclude only rather exotic dark matter models. 

A direct comparison of our results with those of j5j] is 
difficult, because they chose to include in the likelihood 
analysis the substantial statistical uncertainty in the J- 
factor derived from fitting a specific dark-matter profile 
to the stellar-kinematic data. Our limit on (cru) for the, 
e.g., bb channel with DM-mass 100 GeV is a factor 5 bet- 
ter than that of Q , and it is marginally lower than those 
quoted ibidem for Draco, Ursa Major II, Ursa Minor, and 
Goma Berenices. As described in Section HVl using Pass- 
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FIG. 4: 95%-CL upper limits on {av) as a function of the 
neutralino mass for kbb = 1 (curves with triangular markers) 
and k-TT ~ 1 (curves with square markers) for each of the 
four models. The horizontal dashed line shows the canoni- 
cal value of the self-annihilation cross section for a thermal 
WIMP {{av)^ ~ 3 X 10"^^ cm^ s'^). 



7 event selection and the extended source profile would 
increase the upper limit by ~ 40%, whereas the longer 
observing time should decrease it by a similar percentage, 
and thus the effects should roughly balance each other. 

The marginal improvement in the upper limits pre- 
sented here compared with those in ^ quoted for Ursa 
Major II and Coma Berenices is therefore very reason- 
able. However, our nominal J-factor for Segue 1 is a fac- 
tor 2.5 lower than that used in |5|, whereas their quoted 
uncertainty in J is a factor 3.4 which they include in 
the likelihood analysis. Again, one might expect that 
accounting for the uncertainty in J should roughly com- 
pensate the higher nominal J, but our results suggest 
that we would need to set log J = 18.5 to match their 
upper limit, a value that is 2 cr off their nominal J. 

As we have already mentioned, the channels {bb + cc + 
WW + ZZ) dominate for the majority of SUSY models, 
and the total spectrum differs little from that of bb an- 
nihilation. Therefore we may consider the limit given by 
the bb model to be with reasonable accuracy the limit on 
the total cross section. 

But what constraints can the limits place on dark- 
matter models? First, all models for which the dark mat- 
ter has never been in thermodynamic equilibrium should 
predict the annihilation cross section to be much lower 
than the limits shown in Figured If the dark matter was 
once in thermodynamic equilibrium, we can estimate the 
cross section at that epoch from the dark-matter abun- 
dance as {av)a — 3 -10"^^ cm'^ for Majorana particles 
[lll |. As we can see, the upper limits never reach {<Tv)a, 
except for the orbital model; for all other profiles they 
are at least a factor 4 higher. 

On the other hand, the difference is not so large, and 
one might suppose that our results exclude large cross- 
section enhancements, like strong Sommerfeld effect or 



high boost factor C. Unfortunately, even this conclusion 
is, generally speaking, questionable. When dark matter 
was formed, its average speed was tens of thousand of 
kilometers per second. Now it is of the order of hun- 
dreds of kilometers per second [l3l. If the s-channel 
of annihilation dominates (which is typically the case), 
(av) ~ const -I- av^. The p-channel dominates more 
rarely, in which case (cr) ~ const -I- av'^ [20]. We con- 
clude that {(Tv) for small velocities, i.e. for the present 
epoch, can be similar to its value at freeze-out (in the 
first case) or down to 100 times less (in the second case). 
In the latter instance, we cannot even exclude a factor 
of a hundred enhancement. We may only disfavor a very 
large boosting (> 500). 



As we can see, the experimental results obtained do 
not significantly constrain dark-matter models (except 
for the quite weak limitation on the boost-factor or Som- 
merfeld enhancement). A rational question at this point 
is: What observations could improve the limit? First of 
all, a better measurement of the dark-matter profile of 
Segue 1 would reduce the systematic uncertainties (cf. 
Equation [5]), but not necessarily the upper limits. How- 
ever, even then the thermal freeze-out cross-section {av)a 
would not be challenged for any neutralino mass, whereas 
we should reach at least this value to place any reason- 
able constraints on dark-matter models. As we can see 
in Figure 01 the limits on (av) must be improved by two 
orders of magnitude for particle masses around 1 TeV. 
At higher masses the current generation of Cherenkov 
telescopes is more sensitive than Fermi-LAT [2^. For 
satellite-based detectors of GeV-band gamma rays this 
is hardly possible because one would need to operate in 
the background-dominated mode, for which the sensitiv- 
ity increases only with the square root of the observ- 
ing time. The isotropic, presumably extragalactic, back- 
ground alone provides an integrated photon flux (above 
2.4 GeV) per 0.2°-pixel that is 2.5% of our 55-model up- 
per limit. Improvements to the point-spread function of 
the gamma-ray detector would not help on account of 
the intrinsic angular extent of the source. Observations 
of gamma rays with energies ~ 20 GeV may still be con- 
ducted with moderate sky background, although not in 
the photon-counting limit. Therefore, large Cherenkov 
telescopes such as CTA seem to be a promising alterna- 
tive. The main advantage of these instruments is their 
huge effective area ~ 10^ m^ [23], [l^ which could be 
enough to reach {av)a even for heavy WIMPs. The only 
significant weakness of Cherenkov detectors is their rel- 
atively high energy threshold (~ 50 GeV), i.e. they 
are fundamentally unable to detect the annihilation of 
WIMPs of lower mass. On the other hand, for a large 
part of reasonable dark-matter models (like mSUGRA or 
CMSSM) a low-mass WIMP will very soon be excluded 
(or detected) by the LHC. 
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